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Abstract: We describe a new cooling algorithm for SU(2) lattice gauge the- 
ory It has any critical point of the energy or action functional as a fixed 
point. In particular, any number of unstable modes may occur. We also 
^ ' provide insight in the convergence of the cooling algorithms. A number of 

■ solutions will be discussed, in particular the sphalerons for twisted and peri- 
CN . odic boundary conditions which are important for the low-energy dynamics 

of gauge theories. For a unit cubic volume we find a sphaleron energy of resp. 
>• | S s = 34.148(2) and S s = 72.605(2) for the twisted and periodic case. Re- 

markably, the magnetic field for the periodic sphaleron satisfies at all points 
O ' Tr^x = Tr£2 = TiB 2 z . 

m : 

O . 

■ 1 Introduction 

^ ; Saddle points of the energy functional can play an important role in non-perturbative 
dynamics, as they describe the minimal barrier to be taken to go from one to the next 
minimum energy configuration. In a heat bath it determines the critical temperature where 
barrier crossing is no longer suppressed. In a strongly interacting theory it determines the 

. ^ | coupling (or volume) where barrier crossing - quantum mechanically always allowed - is no 
longer exponentially suppressed. The saddle point we have in mind here is the sphaleron 
which has precisely one unstable mode. This comes about since one determines the minimal 
barrier height from a mini-max procedure 0. When one has a path in field space from 
one minimum to the other, one first finds the maximal energy along the path. Then 
one minimizes this maximum over the space of all such paths. The unstable mode will 
correspond to the decrease along the path that goes through this mini-max configuration. 
Given such a saddle point, it is often difficult to prove in all rigour that it is the one with 
lowest energy, as there may be local minima when varying over the space of paths. 

Quite often a sphaleron lies on a particular instanton path (i.e. the action associated 
to the path is also minimal), but this is by no means guaranteed. When the instanton has 
non-trivial parameters (moduli), mini-max with respect to those parameters provides a 
candidate sphaleron. For S 3 one can analytically verify that this candidate is a saddle point 
of the energy functional with one unstable mode |§ . For T 3 this was studied numerically @, 
|[|. It was not entirely conclusive from that study if the sphaleron was indeed associated 
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to the top of an instanton path. With the new algorithm we have been able to address 
this issue more directly and conclude that it is the case for both periodic and twisted 
boundary conditions. 

We will only consider the case of pure SU(2) gauge theories. In an infinite volume the 
mini-max procedure is easily seen to lead to zero sphaleron energy, as instantons can have 
an arbitrarily large scale parameter. In a finite volume, the scale parameter can not be 
larger than the volume, which stabilizes the mini-max procedure at a non-zero energy. In 
the presence of a Higgs field, like in the standard model [0J, the scale is set by the Higgs 
mass. 

Before discussing the numerical results we will first describe the algorithm for finding 
saddle points and discuss the convergence. This study is specific for SU(2), but the principle 
of taking the square of the equations of motion as the action to be minimized will work 
equally well for SU(N) gauge theories, with or without scalar fields. In section 4 we will 
discuss results for the class of analytically known solutions with constant field strength, 
to test the various aspects of the method. In section 5 the results for the sphalerons 
are presented in detail. We end the paper with some concluding remarks. Readers only 
interested in the results should directly go to section 5, perhaps after section 2. 



2 The algorithm 

Consider the square of the equations of motion summed over its variables as a functional. 
Its absolute minima are by construction solutions of the equations of motion. Stability is 
now obvious, but also follows from the fact that at such a solution the Hessian of the new 
functional is (twice) the square of the Hessian of the energy (or action) functional. In the 
n-dimensional continuum we have the following expressions for the two functionals 

8 = -| J d n x Tr(F 2 v (x)) , (1) 

S=-2j d n x Tr((V^ u (x)) 2 ) , (2) 

with F^ u (x) = d^A u (x) — d u A^(x) + [A^x), A u (x)] the Yang-Mills field strength in terms of 
the anti-hermitian vector potential A^x) and T>^ = <9 M + adA^x) the covariant derivative. 

This was considered in the past for constructing monopole solutions in three dimen- 
sions [0], but the process of finding the minimum of S was not gauge invariant. Like for 
the Wilson action || 

S = iV 4 - £ Tr (l - iQ) = N'~ n £ Tr (l - U^x)U v {x + fil/frx + *)Ul(xj) , (3) 

x,n,v Oujl x,fi,v 

a gauge invariant lattice action for the square of the equations of motion is not too difficult 
to write down |J 

S = £ Re Tr(U,(x)Ul(x) - (U,(x)Ul(x)) 2 ) , (4) 

X,fl 

where U^ix) is a group element on the link that runs from x to x + jl (/t is the unit vector 
in the \i direction), iV is the number of lattice points in one of the directions and U^{x) is 
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the sum over 2{n — 1) staples 



= £ + = ^(^(x)C/ M (x+i))C/t(z+/i)+^( x -^)^( x -j));7 !/ (x+/i-i))), 

(5) 

To derive the lattice equations of motion, we observe that 5* depends on Up(x) through 
the expression: 

S(U^x)) = 2N 4 - n ReTi(l-U^x)Ul(x)) . (6) 
As Ufj.(x) is independent of U^(x), it follows that 

U,(x)Ul(x) - U,(x)Ul(x) = . (7) 

Its square, summed over the lattice, yields S. For SU(2) one can also use the fact that the 
equations of motion are solved by 



U fl (x) = ±U fl (x)/\\U fl (x) 



where only the positive sign will allow for solutions that have a smooth continuum limit. 
We can now introduce an error functional Hi 



S = iV 6 -^Tr(l - U,{x)Ul{x)/\\U,{x) 



(9) 



X,fl 



In the continuum limit one finds S = jg(n— 1)~ 2 S. This follows in the usual way [[| from 
expanding the links in a = 1/N, defining the relation between the vector potential in the 
continuum and the link variables on the lattice by 



Uu(x) = Pexp( / AJx + sp)ds) 
Jo 



One finds after some algebra 

U,(x)m(x)/\\U,(x)\\=exp 



2(n - l)N 3 ^ 



YtVuF^x) + 0(N~ 5 ) 



(10) 



We can instead of the standard Wilson action use any improved action, containing 
nxm Wilson loops with appropriate couplings, such as the over-improved actions we have 
used in ref. M 



48 



X,IM,U 



'121 



All that needs to be modified in the above is the definition of which will now also 

contain the staples associated with the larger Wilson loops appearing in eq. (|T2j) 

x - +Mu . Hn x 



U,{x) 



e- 1 
48 



jJL X 



(13) 
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Like for S(e), the e dependence drops out in the continuum limit for S(e). However, for S 
one finds in this limit S(e) = 36((n - 1)(31 - 7e))- 2 S. 

The ordinary cooling algorithm [|H]] for both of these actions is based on the observation 
that the replacement 

U^x) = %(x)/\\U,(x)\\ , (14) 

minimizes S^U^x)). Note that for S (resp. S(e)) cooling we should require the lattice to 
be at least two (resp. three) sites in each direction. To completely specify the algorithm one 
also should prescribe in what order every link is being updated. One sweep is the process 
of updating each link precisely once in that particular order. This defines a mapping 
(U = {U,(x)}) 

U' = T(U) , T = T N oT^ 1 o-..oT 2 oT 1 , (15) 

where Tj is the single-link update of eq. flU] ) and the label i stands for the order in which 
the M = nN n different links are being updated. It is actually not difficult to compute 
exactly by which amount the action is decreasing. For the single-link update we find 

S(Ul(x)) - S(U,(x)) = -2iV 4 -1^(*)|| \\%(x) - U,(x)\\ 2 , (16) 

where we made use of eq. (|I4]). To find the change of the action after one sweep, one has 
to simply add the contributions of each single-link update keeping in mind, however, that 
the value of ||[/ M (a:)|| depends on which links had been updated before. In the continuum 
limit HC/^x)!! — > 2{n — 1) and to a good accuracy the action changes from one sweep to 
the next by the amount 

S(T(U))-S(U) ~ -4(n-l)iV 4 - n ||T(lO-[/|| 2 = -4(n - l)N 4 ~ n ]T ||C^(x) - U^x) || 2 . 

X,fl 

(17) 

For saddle points, which are fixed points of this algorithm, one can of course lower the 
action further and the algorithm is necessarily unstable. 

For the functional S any solution corresponds to the absolute minimum S = 0, and 
cooling with this functional should hence not lead to any instability for the saddle points of 
the original action, eqs. (0,0). Unfortunately, it can be proven that one cannot analytically 
minimize S as a function of one of its links. All that is required, however, is that the 
functional is lowered if and only if we are not at a fixed point. We will see, as our intuition 
tells us, that this is a sufficient condition for convergence of the algorithm (necessity should 
be obvious). Let us nevertheless write down the necessary ingredients to minimize the 
single-link dependence of S, which was already considered by Van der Sijs 0. We restrict 
ourselves in this paper to the version of S derived from the Wilson action, as the analysis 
is prohibitively more complicated for improved actions (for S the situation is in either 
case intractable, because of the appearance of 11^(^)11 in the denominator of eq. ([])). The 
analogue of eq. (§) becomes 

S(U,(x)) = iV 6 -"Re Tr ^(x)W M (x) -S^MWW*)) 2 ~ (^W(*)) 2 ) , (18) 

where the index a runs over the 2(n — 1) different staples of eq. (|5|), such that U^(x) = 
12 a V"(x). The expression for W^(x) is a complicated combination of planar and non- 
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planar Wilson loops, which are best exhibited in terms of unit vectors a, b G {±1, ■ ■ • , ±n} 

w,{x) = 2 y: k 



U a M Q-Jf ^ , 



j*f e m - fn + en - an 



or in explicit form, using the convention that U a (x) = £7_ a (x — a): 

W„(x) = 2Y, a b ^ ±fl . a ^_ b {U a (x)U b (x + a)U„(x + a + S)C/J(a; + a + {i)U\(x + A) 

—U a {x)U fl {x + a)Ub(x + a + p)U^{x + a + b)Ul{x + ajU^x + a)U\(x + //)} 

+ 2 E^^.^W^)^ + + & + /W&V + a + /WlO + A) 

-£/ 6 (x)£/ M (x + S)E/J(a: + /2)£/ a (z + p)U b (x + a + fi)U\{x + b + p)Ul(x + p) 
+Ul(x — a)Ub{x — a)U a {x + b — a)U^{x + b)Ul{x + fi) 

-f/ 6 (x)f/](x + 6-a)^ t (x-a)f/ a (x-a)f/ 6 (x)t/ M (x + S)f/ fe t (x + /i)} . (19) 

The sum in this expression is restricted to those combinations of a and b that do not lead 
to any backtracking in the Wilson loop defined by the products of the links. If one tries 
to minimize S{U^(x)) with respect to U^{x) under the constraint = 1, one finds 

\Up(x) = M(Ur{x))-Wp{x) . (20) 

Note that the lattice should now be at least three sites in each direction. The operator M, 
whose fi and x dependence was suppressed, is defined for SU(2) through the equation 

M(U,(x)) = 6]TTr (u,(x)V«( X y) V«(x) + 2Tr (u,(x)Ul(x)) U,(x) , (21) 

a 

An attempt to solve for A, using the condition that U^(x) is unitary, leads to a complicated 
eighth order polynomial, and solving it numerically for each update would be too costly. 
Instead, we define our update by the simple equation 

r/M _ M(U,(x))-W,(x) 
U ^ X) \\M(U,(x))-W,(x)\\ ' 

The new feature is that the r.h.s. still depends on U^x) and that we do not obtain a 
minimum of eq. flUf). Nevertheless, this is harmless since some algebra will show that our 
update does lower S and as we argued above, and shall show in some more detail in the 
next section, this is sufficient. One finds 



5S(U^x)) = -N 6 - n [\\M(U^x)) - W M (x)|| \\5U fl (x)\\ 2 + < SU^x), M{5U^{x)) >) , (23) 

where, of course, 5U^(x) = U'(x) — U^x), SSQJ^x)) = S(U' (x)) — S^U^x)) and 

<6U,(x),M(6U,(x)) >=3^|Tr(^t (x)V - (x)) | 2 + | Tr( ^t (x) ^ (x)) |2 . (24) 

a 

In the continuum limit one has ||M([/ M (a;)) — VV^x)!! — > 8(n — l)(2n + 1). 
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We note that lattice artefact solutions for S will not be eliminated by eq. (^2|). For 
example, putting one link to minus the identity (while keeping all others equal to the 
identity) is a fixed point of eq. (|22|). It is, however, not a fixed point for eq. (|14D. In 
practise we first use 5* cooling to bring the energy (or action) down to values not too much 
above where we expect a saddle point to occur, after which the S cooling will in general 
bring it to a smooth solution. One can also combine the S and 5* cooling, by replacing 
Ufj,(x) on the r.h.s. of eq. ( p2|) with t/ A1 (x)/||[/ M (x)||. We have not proven that this always 
lowers 5* (although it did in all cases we considered). 

3 Convergence 

We will study the convergence of the algorithm by considering the linear approximation 
for the mapping T in eq. ([15]). For S (or S(e)) cooling, Tj will correspond to eq. (0), 
while for S cooling it will correspond to eq. (^). To warm up we first consider the single- 
link convergence of the latter, while keeping all links (except U^x)) fixed. For notational 
convenience we drop the fi and x dependence. Let us assume that U is a fixed point, and 
X is a small element of the SU(2) Lie-algebra. We now consider U = e x U D and compute 
X' to linear order in X from U' = e x U Q . After some algebra one finds 

X' =\\M(U )-W\\- 1 \6^Tr(XYZ)Y a + 2Tr(XYi)Yj , (25) 

where Y = ^2 a Y a and Y a = Im (U%V a ) = j-(U%V a — V^U ). In the continuum limit this 
is easily seen to imply that ||X'||/||X|| = 0(N~ A ) and we have verified that the single-link 
iteration is indeed extremely efficient. One iteration typically brings U already very close 
to the minimum of eq. ([T5D (although we would like to emphasize again that this is not 
really important, as the other links that contribute to S will be changed during the sweep). 

With respect to the SU(2) Lie-algebra basis ir k /2 {r k the Pauli matrices) eq. fl2"5| ) is 
symmetric, but this is in general no longer true for the linear approximation of T, which we 
will denote by dT. The latter is determined exactly as above by evaluating U' = T(e x U ) 
and extracting X'Jx) to linear order in X from exp(X'(x)) = U'(x)U°{x)\ where U is 
the fixed point. Numerically this is easily achieved by first obtaining the fixed point to 
a high accuracy, assuming the last value of U can be equated to U a , and then looking 
at the response of one complete sweep to U — e x U Q , if necessary doing extrapolation in 
X — > to extract the part linear in X. One might be tempted to prove that T is a so- 
called contraction, for which it is required that ||T(V) — T(i7)|| < k\\V — U\\, with k < 1. 
The contraction theorem asserts that in such a case there is a unique fixed point of the 
mapping T. At the linear level, near the fixed point, this would imply that dT l dT (which 
is symmetric and positive definite) has no eigenvalues bigger or equal to one. We know, 
however, that the fixed point is neither globally nor locally unique. Indeed, for all cases we 
studied there are plenty of eigenvalues bigger than one. If dT would have been symmetric, 
there would have also been eigenvalues of dT that are bigger (in absolute value) than one, 
which would clearly be in contradiction with the fact that the action is lowered with each 
update. 

As dT is not a symmetric matrix, we can not diagonalize it. Still, it can over the complex 
numbers be brought to the so-called Jordan-normal form by a non-singular (complex) 
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transformation A 



dT = A' 1 J A = A' 1 



\ 



A 



(26) 



Ji 



where Jj is the smallest Jordan block, of the form 

(\i i e\ 



Ji 



1 



(27) 



with Aj the eigenvalues. When complex, they occur in conjugate blocks. Different blocks in 
our conventions may have the same Aj, and blocks of dimension 1 correspond to eigenvalues 
with equal left and right eigenvectors. As dT k = A~ 1 J k A, any block with |A$| < 1 will 
approach zero. Blocks with |Aj| > 1 would grow without bound and can never be consistent 
with the fact S is a decreasing functional as a function of the sweeps. What remains are 
the eigenvalues |Aj| = 1. They can, and will, occur for the following reason. Suppose X is 
a zero mode for the Hessian of S at the fixed point U Q , then S(e x U D ) — S(U ) = 0(X 4 ). 
On the other hand using eq. S(e x 'U ) - S(e x Uo) = 0(\\X' - X\\ 2 ), which can only 
be consistent if X' = X. Consequently X is an eigenvector of dT with eigenvalue 1. The 
converse is also true, i.e. an eigenvector with eigenvalue 1 is a zero-mode of the Hessian. 
Consequently, there are at least as many unit eigenvalues as there are gauge degrees of 
freedom. To be precise, there are 3N n minus the dimension of Stab(C/o) unit eigenvalues. 
(Stab(t/o) is the stabilizer of U , which is the subgroup of the gauge group that leaves U 
unchanged. It will only for very special configurations be non-trivial). It is slightly more 
subtle to argue that |Aj| = 1 is only compatible with eq. ( p2]) if Aj = 1 (one uses the fact 
that there can be no linear relation among the eigenvectors, because A is non-singular, 
which implies that after each sweep the amount by which S is lowered is bounded away 
from zero if Aj ^ 1). It is more easy to rule out that Aj = 1 can occur in blocks Jj of 
dimension greater than 1. In other words, left and right eigenvectors for the eigenvalue 
equal to 1 always will coincide. 

Clearly, the fixed point is only unique up to the gauge degrees of freedom and the 
moduli of the fixed point (by moduli we mean the gauge invariant parameter space along 
which S = and consequently along which S is degenerate), and there is a one to one 
relation with the Aj = 1. These do not influence the convergence of the algorithm. But 
it is important to realize that one has no control over which fixed point will be obtained. 
(The gauge in which the fixed point is obtained will entirely be determined by the starting 
configuration, as T commutes with applying a gauge transformation, i.e. T( 9 U) = 9 T(U).) 
It can happen that there are zero-modes of the Hessian that are not associated to a moduli 
parameter of the fixed point, i.e. S is no longer degenerate to fourth (or higher) order. 
In that case one needs to go beyond the linear approximation of T and convergence is 
worse than exponential. In the generic case, however, the rate in which S decreases is 
determined by the gap fi in |A» — 1|, or after k sweeps S decreases the next sweep by an 
amount 5S ~ exp(— 2kfi). The same argument holds of course for S cooling to a (local) 
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minimum, replacing the gap ft by its appropriate value (p) obtained from the mapping T 
for S cooling. 

We can find a relation between the values of these gaps and the minimal non-trivial 
curvature of S at the fixed point if we assume that the gap is due to a real eigenvalue and 
its eigenvector coincides with the fluctuation along the direction of minimal curvature. The 
associated zero-point frequency will be denoted by u (scaled so as to have the appropriate 
continuum limit). At a fixed point U one has for X^ix) = i5A a ^(x)T a / '(2N) 

S(e x U ) = S(U ) + lN- n £ 6Al(x)H£(x,y)6JHXy)+0(6A*) , (28) 

xfia;yub 

which provides the definition for the Hessian H?^(x,y) as a 3nN n square matrix, whose 
eigenvalues will scale to the appropriate continuum limit for N — > oo. Hence, uj 2 equals 
the lowest non-zero eigenvalue (A m ) of the Hessian. It can be easily shown that for any 
saddle point of S, such that S(U ) = 0, the functional S has the expansion 

S(e x U ) = §AT n J2 SAl(x)H; b l/ (x,y)6A b u (y) + 0(5A 3 ) , H = 2H 2 , (29) 

xfia;yub 

where by H 2 we mean of course the matrix multiplication of the Hessian of S with itself. 
In the same shorthand notation (i.e. dropping the indices), we assume that under the S 
cooling H ■ X = X m X and X' = dT ■ X = (1 — p)X, which allows us to calculate the 
change in the action both from eq. (fZS|), giving 5S = —4u 2 p(l — ±/i)iV 2 ~ n ||X|| 2 and from 
eq. flTTD , giving SS = —A{n — l)/i 2 iy 4_ri ||X|| 2 . Equating these results gives 

M ~ (n - l)iV 2 ~ (n - 1)A^ 2 ' 1 ' 



Similarly, using eqs. (|23|J29| ) one finds for cooling with the functional S 

(n-l)(2n + l)N 4 ~ (n - l)(2n + 1)N 4 ' 1 J 

This is of course in accordance with ones intuition that the algorithm can slow down 
considerably if the functional has directions in which it is very shallow. Note that our 
results indicate that slowing down of the algorithm goes proportional to N 2+n for S cooling 
and proportional to N 4+n for S cooling (in both cases a factor N n comes from the number 
of links to be updated in a single sweep). 

Finally, most of the moduli for a continuum solution are removed by lattice artefacts, 
although we will discuss in the next section that in some cases a few of them can survive 
on a lattice. An important example of moduli that are removed by lattice artefacts are 
those that are related to the translational invariance of the continuum theory on a torus. 
As this breaking is due to the fact that integrals in the continuum are replaced by sums on 
the lattice, one can easily argue, as in ref. ||, that the breaking is very mild and goes as 
exp(— pN) where p is the typical scale of the continuum solution in units of the length L of 
the torus. For small lattices these translational modes can cause some trouble, but those 
lattices are usually anyhow too small to extract reliable information about the continuum 
limit. In ref. [ffl we made judicious use of the dependence of the lattice action on the 
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moduli parameters, by using the over- improved action S(e = — 1) (see eq. (0)) so as to 
ensure that the instanton with the maximal scale parameter allowed by the finite volume 
is obtained as a local minimum for the action. With S cooling we could then find rather 
accurate lattice solutions, whose scaling to a continuum solution was studied in detail. 
One of the advantages of the S cooling is that one can verify if the fixed point has been 
reached, as S should approach zero. This can, of course, also be used as a check under 
S cooling. For example, we performed for an instanton on an 8 3 x 24 lattice cooling for 
about 6600 sweeps, which allowed for a very precise exponential fit to the decay of the 
action. The effective value for \i deduced from this fit is 9.8 x 10~ 5 (see fig. 3 of ref. ||). 
We find for the best configuration of this four dimensional instanton (see table I of ref. [f|) 
S(e = -1) = 4.8 x 10" 4 and for S(e = -1) = 1.3 x 10~ 6 . The results obtained with 
S(e = 0) cooling on this configuration (see the same table), lead to S(e = 0) = 2.0 x 10~ 6 , 
S(e = 0) = 8.2 x 10~ 9 , and effectively fi ~ 10~ 6 (which can be understood since for S(e = 0) 
the 0(a 2 ) corrections to the action are canceled, removing to this order the dependence of 
the action on the scale parameter. Using the heuristic estimates of \i discussed above we 
find fj, ~ u 2 /(3N 2 ). Since u(e ^ 0) = 0(N~ l ), whereas u{e = 0) = 0(N~ 2 ), one obtains 
with N = 8 the right orders of magnitude for fi). In the remainder we will discuss the new 
results for n = 3, i.e. results for the Yang-Mills energy functional. 



4 Results: Constant solutions 

We will first describe the case of constant curvature solutions, which are known analytically 
both on the lattice and in the continuum. We take from ref. [|11[| the following results. If 
we introduce the so-called twist tensor n^, which is antisymmetric and takes integer values 
(equivalently specified by the magnetic flux = ^€ij k rij k ), the following configuration is a 
solution on a symmetric three torus of length L — 1, or a symmetric lattice with N 3 sites 



A°Ax) = i(-irn jk x k + C,-)? , U?(x) = exp(A°( a ;/iV)/iV) (32) 

2 



These fields are periodic up to a gauge transformation, to be specific [|12|] 

Aj(x + k) = n k (x)(A J (x)+d j )Q^ 1 (x) , Uj(x + Nk) = n k (x/N)U 3 (x)n k -\(x + ])/N) , 

(33) 

where 

Q k {x) = exp^inxjnjkTs) . (34) 

If some of the components of the twist matrix are odd, quantities that transform in the 
fundamental representation of SU(2) are anti-periodic, which can be most easily seen from 
the Polyakov line 

Pj(x) = |Tr{Pexp(^ ds A^x + sj))Qj(x)} = cosQQ - im jk x k ) . (35) 

The identical result follows on the lattice, after replacing x by x/N. For the continuum 
Yang-Mills and lattice Wilson actions (energies) we find 

S = - 2 E4 , S = 4iV 4 £sin 2 g|) , (36) 
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which are constant and independent of C. The latter form the gauge invariant moduli of 
these solutions, both on a lattice and in the continuum. If the action is non-zero, i.e. the 
twist is non-trivial, two of the moduli can also be identified with translational invariance, 
whereas translation in the other independent direction (which forms a zero eigenvector of 
the twist matrix) leaves the configuration invariant (on the lattice there is of course no 
continuous translation invariance). The stabilizer of this configuration is the subgroup 
of constant abelian gauge transformations, which is one dimensional. On the lattice this 
means that the Hessian has 3N 3 + 2 exact zero-modes and that dT has the same number of 
unit eigenvalues (which we verified numerically). In the continuum we can even compute 
the Hessian exactly ||11|| . The eigen modes are products of a theta function (compare the 
study in the context of the "Copenhagen vacuum" |Tj|) and a plane wave (purely plane 



waves for abelian fluctuations). One easily shows that the gauge field is (anti) periodic in 
the fh direction, i.e. A(x + fh/e) = t^A(x)t^ a , with e = |g.c.d.(mj)| and A(m) = or 1, 
its value being determined from (repeated use of) eq. (0). The spectrum is now given by 

A±(n,f) = 27r|H|(2n + l±2)+e 2 ((2£ + A)tt + C ■ m/e) 2 /\\m\\ 2 , X (k) = (2tt£) 2 , (37) 

The multiplicities are 2|e| for X± and 2 for A D . The parameters (n, i ,k) run over 
(IN, 7L, 7L Z ). On the lattice \ {k) will be replaced by AN 2 J2i sin 2 (nki/N), but the effect of 
the lattice on A± is not expected to admit an analytic expression. In fig. la-c we plot the 
first few eigenvalues of the Hessian for resp. fh = (0, 0, 1), fh = (0, 0, 2) and fh = (1, 1, 1), 
computed numerically for N = 4 and 6 in comparison with the exact continuum result as 
a function of C, of which only the component along fh is relevant. Only for fh = (1, 1, 1) 
is A non-trivial. 

A new feature that we believe not to have been observed so far, is that there are constant 
(^ 0) curvature solutions on a torus that are marginally stable. To be precise, for ||m|| = 1 
and \C ■ fh — n\ < ir — \^2n (using the 2n periodicity of the spectrum in C ■ fh we can 
restrict to C ■ fh G [0,27r]) there are no negative eigenvalues for the Hessian. We call this 
marginal stability, as one can first bring C ■ fh outside the above range without changing 
the energy, and then let the configuration decay along one of the two unstable modes. For 
||m|| 7^ 1 one can easily see that this effect doesn't occur. 

Nevertheless, in all cases considered in fig. 1 the Hessian can have accidental zero 
eigenvalues for particular values of C, such that potentially the algorithm can slow down 
dramatically. In table I we compare the gap fi (the distance of the largest non-unit eigen- 
value of dT to 1) with the predicted value fi = A 2 n /(14A^ 4 ), determined from the smallest 
non-zero (absolute) eigenvalue A m of the Hessian for S (see eq. (|3l|) . The agreement is (per- 



haps surprisingly) very good. For fh = (1, 1, 1) and A^ = 4 there is considerable slowing 
down. What is called slowing down is of course dependent on the available computer power. 
We consider something slowed down considerably if (several) 10.000 or more sweeps are re- 
quired to half the value of S. In the case of slowing down, we could generally achieve good 
results by performing a small number of S coolings, when the S cooling is being "stuck", 
sometimes after a (very) small random update on the links. That S cooling can sometimes 
help can be argued for the case that A m > as follows. We have seen that to some degree 
the A m eigenvector is also an eigenvector of dT, corresponding to the eigenvalue A = 1 — fi. 
This should be the dominant mode along which the configuration differs from the true fixed 
point. All other modes, in particular the unstable ones, should have decreased significantly 
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under the S cooling. If we don't cool too long with S, the unstable mode does not get a 
chance to grow, whereas the direction of the stable mode keeps on contracting (faster by 
a factor 14iV 2 /A m ). With S cooling one can re-contract the unstable mode, which might 
have been grown during the S cooling. We expect that the extra computational effort will 
not exceed by much what one would have to spend in trying to accelerate the algorithm, 
e.g. using fast-fourier [[14], [7|] or multi-grid |j| methods. 



N 


fh 


C ■ fh 


S 


S 


5th 


A m 




Ath 


4 
6 
8 


(1,1,1) 


3.472 
2.517 
3.175 


2.90 ■ 10" 4 
1.60 • 10- 2 
2.65 • lO" 6 


59.02774641 
59.17251733 
59.20573728 


59.02761860 
59.18005528 
59.20573660 


0.051 
-0.574 
1.906 


6.49 ■ 10~ 7 
1.74- 10~ 5 


7.31 ■ 10~ 7 
1.82 • 10~ 5 
6.34- 10~ 5 


6 


(0,0,2) 


4.229 


5.24- 10~ 7 


78.75660845 


78.75660838 


4.084 


8.58 • 10~ 4 


9.19 • 10~ 4 



Table I: Results for the constant curvature solution on a iV 3 lattice with 
magnetic flux fh. Sth is the value obtained from the exact lattice result given 
in eq. flSTf), whereas /} t h = ^ m /(14-/V 4 ), the predicted value for the distance 
of the largest non-unit eigenvalue of dT to unity. The smallest (in absolute 
value) non-zero eigenvalue for the Hessian of 5* is indicated by X m . This, as 
well as /t, was obtained from the exact lattice solution with the tabulated 
values of C ■ fh extracted from the data. 



Not only 5* itself is a good indication for how far one is from the fixed point. Also 
dT can, and usually does, have an eigenvalue bigger than 1 if one is not close enough to 
the fixed point (the gauge modes always correspond to unit eigenvalues), whereas only 
at a fixed point will the gauge modes be exact zero-modes for the Hessian. All this 
provides us with more than enough criteria to be sure that we have an accurate lattice 
solution of the equations of motion. In table I we also compare the energy of the constant 
curvature solutions with the exact lattice result of eq. fl3"Tj|). The agreement is excellent, 
except perhaps for N = 6 and m = (1, 1, 1), but in that case we did not try to cool the 
configuration down to a very small value of S. Nevertheless, the error in the energy is 
smaller than one part in 6000. Also quantities like the Polyakov loops agree very well with 
the exact lattice expressions, eq. (135|) . 



5 Results: The sphaleron 

We now discuss the results for the new sphaleron solutions. Some of the results are collected 
in table II. Also here we see that convergence of the 5* cooling algorithm can be well 
predicted in terms of A m . In figs. 2a-c we show scatter plots for the (complex) eigenvalues 
of dT for the fh = sphaleron, using a checkerboard updating for fig. 2a (to allow for 
vectorization of the algorithm) and a sequential updating for fig. 2b, both at N = 4. We 
see that the checkerboard updating makes the bulk of the modes contract much faster (in 
addition to the gain of computational speed due to the vectorization). This remains true 
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on larger lattices. In fig. 2c we show for comparison the checkerboard result for N = 8. The 
eigenvalues that are responsible for the slowing down of the algorithm near the fixed point, 
in particular the eigenvalue closest to 1, will only very weakly depend on the order in which 
one updates the links. Table II shows that the sphaleron with periodic boundary conditions 
(m = 0) where X m ~ 1.3, has considerable slower cooling than the sphaleron with twisted 
boundary conditions, for which X m ~ —10.8. Nevertheless, we managed to obtain rather 
accurate results. In table II we also list S, S and the blocked Wilson action S2X2, for which 
the lattice spacing is effectively twice as large. Like the instantons of maximal size studied 
in ref. ||, [|, the sphaleron is a very smooth solution. We do not require huge lattices to 
approach the continuum limit. Together with easily obtainable expansions in the lattice 
spacing a = 1/N, we have a large enough window (available values for N) to get accurate 
results. We can extract, as for the instantons (see ref. Q for details), the continuum action 
by fitting to the formula 

Snxn = S (l- an 2 /N 2 - (/V + in 2 + 5)/N 4 ) + 0(n 6 /N 6 ) . (38) 

We find for the twisted continuum sphaleron (fh = (1, 1, 1)) an energy £ = 34.148(2) and 
for the periodic one (m = 0) we find £ = 72.605(2), with conservative error estimates. 



N 


in 


Am 






s 


5'2x2 


S 


4 


(1,1,1) 


-11.128 


2.73 


■ 10" 


-2 


3.46 


■ io- 


-2 


33.24878 


29.77153 


7.90 


■ io- 


-6 


6 




-10.900 


8.74 


• io- 


-3 


6.55 


• io- 


-3 


33.77789 


32.49357 


2.28 


• io- 


-7 


8 




-10.817 








2.04 


■ io- 


-3 


33.94548 


33.28093 


1.20 


• io- 


-11 


12 












4.03 


• io- 


-4 


34.05973 


33.78353 


2.34 


• io- 


-6 


4 


(0,0,0) 


1.704 


7.33 


■ io- 


-4 


8.10 


■ io- 


-4 


68.13887 


52.59651 


2.73 


■ io- 


-5 


6 




1.415 


1.08 


• io- 


-4 


1.10 


• 10" 


-4 


70.70174 


64.36296 


2.11 


• io- 


-5 


8 




1.301 


2.41 


■ 10- 


-5 


2.95 


■ io- 


-5 


71.55060 


68.16790 


6.93 


■ io- 


-9 


12 












5.83 


■ io- 


-6 


72.14095 


70.70631 


4.13 


■ io- 


-5 



Table II: Results for the sphalerons on a iV 3 lattice with twisted and periodic 
boundary conditions. The smallest (in absolute value) non-zero eigenvalue 
for the Hessian of S is indicated by \ m . The predicted value for the distance 
of the largest non-unit eigenvalue of dT to unity is given by £i t h = A^/ (14iV 4 ). 
We also give the value of the Wilson action, its blocked value ^2x2 and S. 



The evidence that these configurations are sphalerons follows from a study of the eigen- 
values for the Hessian of S, whose low-lying values are listed in table III for iV = 4, 6 and 
8. We have verified that there are exactly 3N 3 + 3 zero- modes for both sphalerons, corre- 
sponding to the gauge modes and the three translational moduli (for N > 6 the latter are 
indistinguishable from the gauge zero-modes, for N = 4 their eigenvalues are not bigger 
than 0.01). Most importantly we have verified that there is precisely one unstable mode. 
Furthermore, we have not been able to find any other saddle point with a lower (non-zero) 
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energy than the energies of the present configurations. For all practical purposes we thus 
consider the newly found solutions to be sphalerons. 

We now wish to show that like on S 3 || , these sphalerons are at the top of a tunnelling 
path. For the case of twisted boundary conditions this is the easiest to demonstrate, since 
the instanton with a continuum action of 4ir 2 , which was extensively studied on the lattice 
in ref. ||, has only the trivial translational and some discrete moduli. This means we 
only have to demonstrate that the top of the instanton path is identical to the sphaleron. 
There are two reasons why the top of an instanton path might not correspond exactly to 
a sphaleron. First, the lattice has a finite resolution and second, the extent in the time 
direction is finite. In particular the finite resolution causes problems because one needs to 
extrapolate along the tunnelling path the value of the energy and of S. In particular the 
latter has a sharp dip near the top. A very large value of A" would be required to obtain a 
sufficient resolution. The effect of the finite extent in time is expected to be small since the 
instanton path approaches the vacua exponentially in time 0. In fig. 3 we plot S ~ S/64 
for each time-slice along the tunnelling path, cmp. ref. ||. 



m = (1, 1, 1) sphaleron m = sphaleron index 2 



N — A 


N 

1 V 




N — R 

1 V — o 




N — A 

1 V l ± 


N — (\ 

1 V — u 


N — R 

1 V — O 




N 

1 V 


— Q 

— o 


-11.1280 


-10 


9003 


-10.8174 




-13.1960 


-13.3935 


-13.4432 




-13 


1170 


0.0074 





0001 


0.0000 




0.0073 


0.0006 


0.0000 




- 


7466 


0.0091 





0001 


0.0000 




0.0124 


0.0010 


0.0000 







0000 


0.0101 





0001 


0.0000 




0.0139 


0.0011 


0.0000 







0000 


13.1667 


13 


3371 


13.3934 




1.7041 


1.4151 


1.3009 







0000 


13.1667 


13 


3371 


13.3934 




1.7070 


1.4303 


1.3012 




3 


2618 


13.1667 


13 


3372 


13.3934 




7.9677 


8.0984 


8.1260 




5 


8733 


16.0475 


16 


8663 


17.1597 




7.9682 


8.1008 


8.1261 




5 


8807 


16.0478 


16 


8663 


17.1597 




7.9684 


8.1279 


8.1266 




6 


1738 


16.0484 


16 


8663 


17.1597 




8.9984 


10.2977 


10.6604 




8 


5986 


21.4686 


22 


9444 


23.4481 




9.0189 


10.3306 


10.6611 




8 


6069 


21.4747 


22 


9444 


23.4481 




9.9919 


10.8287 


11.0701 




15 


0281 


21.8119 


23 


6594 


24.2198 




10.0256 


10.8312 


11.0701 




16 


7321 


21.8260 


23 


6594 


24.2198 




10.0354 


10.8594 


11.0707 




19 


7489 


25.0998 


26 


1760 


26.4962 




24.6675 


26.4047 


26.9618 




20 


6039 


25.1050 


26 


1760 


26.4962 




24.6703 


26.4205 


26.9621 




21 


3605 



Table III: The lowest lying eigenvalues for the Hessian of S for the twisted 
and periodic sphaleron on a A^ 3 lattice for N = A, 6 and 8. Also is given the 
spectrum for the periodic solution with two unstable modes (index 2), close 
in energy to the sphaleron, discussed in the text (see also table IV). 
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We take the configuration at the top of the instanton path, where S is minimal, and cool 
it with S. This way we easily obtain the sphaleron solution. As S decreases monotonically 
under the S cooling, and as S is already small at the top of the instanton path, we 
expect that the instanton path goes through a sphaleron. Further evidence is provided by 
calculating the quantity 

D N (t) = qu t - U s \\ 2 /N ~ J d 3 x(5A«(x)) 2 , (39) 

where the latter expression is obtained in the continuum limit. Ut is the configuration at the 
time t along the instanton path and U s is the sphaleron configuration, obtained after cooling 
from the configuration U t (for this, t should be chosen not too far away from the top of the 
barrier, as otherwise cooling will bring the configuration to the vacuum). As the cooling 
algorithm commutes with gauge transformations, this gives a gauge invariant measure for 
the distance of the two configurations. Indeed, computing \\U t — U t >\\ gives huge values, as 
the gauge between different time slices will differ randomly. For the configurations of fig. 3 
we find D B (0) = 0.0246, D 8 (l/8) = 1.162, D g (-l/8) = 1.5713 and D 12 (0) = 0.0340. We 
have indicated also the values at either side of the top of the barrier for N = 8 to show 
that, not being exactly at the top can have a sizable effect on -Djv(O) (and on S(t = 0)). 
We believe the numerical evidence is sufficiently convincing to conclude that the instanton 
path does go through a sphaleron. In fig. 4 we give for iV = 12 the two dimensional cross 
sections of £b(x,i/,z) at the twelve lattice values of x. Because of the cubic symmetry of 
the sphaleron, the other cross-sections have identical profiles (modulo translations). 

For the case of periodic boundary conditions the instantons have additional moduli. 
We did provide evidence in ref. |4| that for T — ► oo the moduli are described by the 
scale, translational and vacuum parameters. The latter being specified by the Polyakov 
loops Pi(t = oo) = —Pi(t = — oo) (we have not yet proven or disproven if this relation 
at t — ±oo can be relaxed). In particular we probed the dependence of the tunnelling 
path on P 3 = P 3 (t = — oo), keeping the other two equal to unity. In practise we can of 
course not reach \t\ = oo, and this is a problem for P% = ±1. In that case one can show 
that the approach to the vacuum for large t is no longer exponential due to the quartic 
behaviour of the potential in the zero-momentum sector. We demonstrate the sensitivity of 
the instanton path on the extent in the time direction for N = 8 and N t = 24, 48 in fig. 5. 
We see that for large N t the tails in the electric energy [|], |5| are lowered considerably, while 
the magnetic energy at the top of the barrier (not visible at the scale of fig. 5) is lowered 
from S(e = -1) = 74.041 at N t = 24 to 73.465 at N t = 48, actually below the value of 
73.923 obtained for P 3 = 0.5 in ref. || (for which we did not observe such a sensitivity 
in N t ). It means that there is near degeneracy in energy at the top of the barrier as a 
function of P 3 . The variation of the energy as a function of P 3 can be so small, because the 
energy functional near the sphaleron has two rather shallow directions (associated to the 
two eigen modes with A m ~ 1.3). To investigate if our sphaleron indeed corresponds to the 
top of a tunnelling path, presumably the one with P% = 1, we took the configuration at the 
top of the energy barrier for the values of P 3 that were considered in ref. |J and applied 
S cooling for 8000 sweeps on each, such that the energy and S stabilized, and further 
cooling would make the configuration only move (very slowly) along the slowest mode in 
the cooling algorithm. We discovered a very useful set of parameters to distinguish the 
various configurations thus obtained. It turns out that the energy profile is almost identical 
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for each configuration, that the profiles for the Polyakov lines only vary weakly, but that 
the vector 

b= (Tr(P 2 ),Tr(P 2 2 ),Tr(P 3 2 )) /^£Tr(P 2 )* , (40) 

averaged over the lattice (< b >), clearly differentiates the various configurations. The 
sphaleron solution distinguishes itself by a symmetric value b = (1, 1, 1) at each point 
independently. In all other cases the variance of the parameter b is non-zero. To a good 
degree < b\ >=< b 2 >, which is of course a direct consequence of the symmetry of the 
lattice and of the fact that P\ — P 2 . We also compute D' 8 , which is defined as in eq. (|3~9"D 
but with Ut the configuration obtained after the S cooling. For this it is important to note 
that the instantons in table IV were all generated from the same random start at P3 = 1 
(or C 3 = 0; see ref. This means that each is in the same gauge and D' 8 is a good gauge 
invariant measure for the distance of the various cooled configurations to the sphaleron. 
In table IV we collect the results thus obtained, and see that after 8000 cooling sweeps S 
is dramatically lowered also for P3 = 0. If we compute the Hessian for this configuration 
we find a Morse index of two, i.e. there are two unstable modes (see table III for the 
low- lying eigenvalues). Note that S for this configuration is only slightly higher than for 
the sphaleron, and that X m = —0.74655 will cause an even slower cooling than for the 
sphaleron. We have performed 8000 extra S cooling sweeps before computing the Hessian, 
but the quality of its zero modes were sufficiently good to trust the cooled configuration 
for P 3 = to be close enough to a solution. We notice that except near P 3 = 1, the value 
of < b > at the top of the barrier (denoted by < b t > in table IV) is changed after 8000 
S cooling sweeps. We conclude that the index 2 solution is not at the top of any of the 
instanton paths we considered here. 
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<b t > 


< b > 


S 


5*2x2 


S 


D's 





(1.01,1.01,0.99) 


(1.01,1.00,0.99) 


71.5507 


6£ 


5.1680 


2.08 


• 10" 


-4 


0.000 


tt/5 


(1.02,1.02,0.96) 


(1.03,1.02,0.95) 


71.5547 


6£ 


5.1732 


9.77 


• 10- 


-3 


0.056 


2tt/5 


(1.05,1.05,0.89) 


(1.07,1.06,0.85) 


71.5913 


6£ 


5.2205 


7.53 


• 10- 


-2 


0.100 


2tt/3 


(1.10,1.10,0.77) 


(1.14,1.13,0.65) 


71.7569 


6£ 


5.4322 


1.58 


• 10- 


-1 


0.560 


4tt/5 


(1.11,1.11,0.73) 


(1.15,1.15,0.60) 


71.8301 


6£ 


5.5250 


1.18 


• 10- 


-1 


0.894 


71 


(1.12,1.12,0.69) 


(1.17,1.16,0.53) 


71.9318 


6£ 


5.6537 


2.60 


• 10- 


-4 


2.186 



Table IV: Results for N = 8 associated to the top of the instanton path 
(obtained with S(e = —1) cooling, see ref. ||) that interpolates from the 
vacuum specified by Pi = P 2 = 1 and P3 = cos(C*3/2) at t — —00, to the 
vacuum specified by P\ = P 2 = —1 and P 3 = — cos(C3/2) at t — 00. The 
vector b is the parameter defined in eq. fl40|). From the third column on are 
tabulated the results after 8000 S cooling sweeps. D' 8 is defined as in eq. (p9f) , 
with Ut the cooled configuration. Both from the first and last case one finds 
solutions of the equations of motion, resp. the sphaleron and the index 2 
saddle point. See table III for their Hessians. 
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Finally, we calculated .Djv(t) from the configuration in table IV with C3 = (P3 = 1) 
after further S(e = — 1) cooling on the instanton, as represented in fig. 5 (at N t = 24). 
We find D 8 (0) = 0.0496, D 8 (-l/8) = 2.212, D 8 (l/8) = 3.369 and D 12 (0) = 0.0368 (not 
represented in fig. 5). Taking into account the time-tail effects discussed before and the 
fact that the instanton path was generated with S(e = —1) cooling, whereas S cooling was 
performed at e — 1, the quality of these results are sufficient to conclude that the sphaleron 
lies on the instanton path that has Pi(t = —00) = —Pi(t = 00) = ±1. In fig. 6 we give the 
results for the energy distribution on a 12 3 lattice. 



6 Conclusions 

The main motivation for this paper was to find the sphaleron solution for pure SU(2) gauge 
theories in a finite volume and to show that it is at the top of the barrier that separates 
two classical minima, connected by an instanton solution with minimal action. We studied 
both the case of twisted and periodic boundary conditions. For the latter the situation 
was rather subtle because of the many additional moduli in the instanton parameter space. 
Furthermore, at the sphaleron there are two very flat directions for the energy functional, 
and we found a saddle point with two unstable modes not more than 5S = 0.5 above the 
sphaleron. The sphaleron solution for the periodic case has some very special properties. 
Like for the twisted case |J only low fourier components seem to dominate, but especially 
the property that locally Tr(Bi) 2 is independent of i is intriguing. There is good hope 
that in the future this problem might be tractable analytically. In any case it presents an 
interesting challenge. 

To achieve all these results we built on the idea that one can find unstable solutions to 
the equations of motion by minimizing the functional obtained by squaring the gradient of 
the field equations J7|, |9|. What we have added is a specific algorithm for SU(2) to minimize 
this functional in a deterministic way. This allowed us to analyse in detail the convergence 
of the algorithm, also providing results for the standard cooling algorithm [ID|. It is likely 



that the algorithm can be extended to include a scalar sector so that we can also study 
the sphaleron in the standard model [0. In any case one can always try to minimize the 
functional by doing random updates. One can easily think of many more applications 
worthwhile pursuing with these methods. 
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Figure captions 



Figure 1: The low- lying eigenvalues of the Hessians for the constant curvature 
solutions specified by (a) m = (0, 0, 1), (b) m = (0, 0, 2) and (c) m = (1, 1, 1) 
at iV = 4 (dotted lines), N = 6 (dashed lines) and N = oo (i.e. the con- 
tinuum). Each level is two-fold degenerate except for rh = (0,0,2) and the 
levels indicated by the arrow in (a), which are all four-fold degenerate. 

Figure 2: The scatter plots for the complex eigenvalues of dT for the periodic 
sphaleron at N = 4 (a,b) and N = 8 (c), for the latter we also show the 
blown- up region near A = 1. Figs. (a,c) were generated with a checkerboard 
updating, whereas fig. (b) was generated with sequential updating. 

Figure 3: The value of Sit) ~ S(t)/64 (at e = 1) for the twisted instanton M 
with continuum action 47r 2 on a 8 3 x 24 (triangles) and a 12 3 x 36 (squares) 
lattice. The strong dip in S(t) at t — indicates stationarity of the energy 
functional due to the tunneling through a sphaleron. 

Figure 4: The profile for the magnetic energies £ B for the twisted sphaleron 
on a 12 3 lattice at fixed values of function of y and z. The lattice 

is represented by the grid-lines (all coordinates are interchangeable). The 
vertical axes runs from to 110, in units of the inverse physical length of the 
box. 



Figure 5: The electric and magnetic energies obtained after cooling with 
S(e = — 1) for a periodic 8 3 lattice in the space directions and for boundary 
conditions fixed in the time direction of N t sites to Pj = — 1 at one end and 
Pj = 1 at the other end, for N t = 24 (squares for £ E and crosses for E B ) 
and N t = 48 (triangles for £ E and stars for £ B ). The time-tails are plotted 
at a blown-up scale to show the electric tails || that distort the instanton 
solution at small Nt- 

Figure 6: The profile for the magnetic energies £ B for the periodic sphaleron 
on a 12 3 lattice at fixed values of function of y and z. The lattice 

is represented by the grid-lines (all coordinates are interchangeable). The 
vertical axes runs from to 110, in units of the inverse physical length of the 
box. 
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